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ABSTRACT 

Reverberation lags have recently been discovered in a handful of nearby, variable AGN. Here, we analyze a 
^100 ksec archival XMM-Newton observation of the highly variable AGN, ESO 1 13-G010 in order to search 
for lags between hard, 1.5 - 4.5 keV, and soft, 0.3 - 0.9 keV, energy X-ray bands. At the lowest frequencies 
available in the lightcurve (< 1 .5 x 10 -4 Hz), we find hard lags where the power-law dominated hard band lags 
the soft band (where the reflection fraction is high). However, at higher frequencies in the range (2-3) x 10~ 4 
Hz we find a soft lag of -325 ± 89 s. The general evolution from hard to soft lags as the frequency increases 
is similar to other AGN where soft lags have been detected. We interpret this soft lag as due to reverberation 
from the accretion disk, with the reflection component responding to variability from the X-ray corona. For a 
black hole mass of 7 x 10 6 M Q this corresponds to a light-crossing time of ^9 R g /c, however, dilution effects 
mean that the intrinsic lag is likely longer than this. Based on recent black hole mass-scaling for lag properties, 
the lag amplitude and frequency are more consistent with a black hole a few times more massive than the best 
estimates, though flux-dependent effects could easily add scatter this large. 

Subject headings: accretion, accretion disks — galaxies: active — galaxies: nuclei — X-rays: galaxies 



1. INTRODUCTION 

Since the discovery of t he first X-ray reve rberation lag 
in the AGN 1H 0707-495 dFabian et alll2009l) . X-ray lags 
between lightcurves from the soft excess and power-law 
dominated regions of the spectrum ('soft ' lags) have now 
been observed in over a dozen AGN dZoehbi & F abian 
201 U lEmmanoulopoulos et alJ I20TH; iTripathi et alj 120111: 



De Marco et alj|201 U 120121; iFabian et al J 120 121) ancfalso i n 
the black hole X-ray binary GX 339-4 dUttlev etalJl20Tfl) . 
Recently, the first disco very of a broad Fe K lag has now been 
reported in NGC 4151 (IZoghbi et al.ll2012h . with large XMM- 
Newton datasets from 1H 707-495 and IRAS 13224-3809 
also revealing an Fe K lag dKara et al.ll2012alfbT) . Such lags 
are naturally predicted in the reflection paradigm where a 
corona (the power-law component) irradiates the accretion 
disk, leading to a reflected component consisting of fluores- 
cent emission lines and scattered continuum emission that are 
broadened and skewed by the dynamical and relativistic ef- 
fects present in the accretion disk close to the black hole (see 
iRevnolds & Nowakll200a lMiller||2007l for reviews of reflec- 
tion). The lags can then be ascribed to the light travel time 
between the dire ct power-law compon ent and the reflection 
components (see IRevnolds et alJll999L for a theoretical dis- 
cussion on X-ray reverberation). 

Lags in the opposite sense ('hard' lags where the soft 
band leads the hard) are also seen on longer timescales. 
These hard lags have been seen in both X-ray binaries and 
AGN (e.g. iMivamoto & Kitamotdll989l: iNowak et al.l ll999fc 
Papadakis et al. 2001; Arevalo et al] l200a) . and are thought 
to arise due to viscous propagation of mas s accretion fluctua- 
tions in the disk transmitted to the corona (iKotov et aUlfOOlt 
I Arevalo & Uttlevll2006t lUttlev et ail 1201 lh . Comparing the 
evolution of the lags with Fourier frequency from all the AGN 
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where they have been detected, a general shape emerges, with 
hard lags at the lowest frequencies, transitioning to soft lags at 
higher frequencies, interpreted as viscous propagations dom- 
inating on long timescales and reverberation dominating on 
shorte r timescales. This i s not the only proposed model for the 
lags -[Miller et al. (2010) suggest an alternative interpretation 
involving reflection from distant absorbers homogeneously 
distributed along our line of sight, something that requires a 
specific geometry. While such a model could potentially ex- 
plain one source, the presence of soft lags in a large number 
of AGN rules out this model in general (further arguments 
against this model are given in Zoghbi et al. 20TJ). 

In this manuscript, we present an X-ray timing anal- 
ysis of ESO 113-G010 (z = 0.0257), using XMM- 
Afewfonobservations, which we find to show the same lag- 
frequency evolution as has been observed in other AGN. 
ESO 113-G010 is a Seyfe rt 1 .8 displaying a strong soft ex- 
cess and X- ray variability dPietsch et al.lll998c iPorquet et alj 
12004 120071) which are the two basic components needed for 
tracing negative reverberation time delays. 

2. DATA ANALYSIS 

ESO 113-G010 was observed by XMM-Newton on 
2005/1 1/10 for a to tal of 104 ksec, ObsID: 0301890101 (see 
IPorquet eTa l. 2007, for a previous analysis of these data). We 
analyzed the XMM-Newton data using S AS version 1 1 .0.0 and 
the most recent calibration files. Calibrated events files were 
generated from the Observation Data Files using the epproc 
and emproc tools for the EPIC/pn and EPIC/MOS detectors, 
respectively. We searched for background flares in the usual 
manner, looking at the lightcurves from the entire detector 
above 10 keV for the MOS, and in the range 10-12 keV 
for the pn. We found a significantly elevated background 
count rate during the beginning and end of the observation, 
which after we filtered out yielded a continuously sampled 
data stream (important for timing analysis) of 92ks. 
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FIG. 1.— XMM-Newton/pn spectrum of ESO 113-G010 fitted with a 
power-law (blue, dashed line), blurred reflection (dotted, red line) and two 
narrow Gaussian emission lines (green, dot-dashed line). 



TABLE 1 
Spectral fits to ESO 113-G010 



Parameter 


Value 




N H (i0 20 cm" 2 ) 


2.74 (fixed) 




^cdgc (keV) 


0.86 ±0.01 






0.26 ±0.01 




Eu„, CkeVI 


6 52 + £° 3 




fline (keV) 


0.0 (fixed) 




Nline 


(1.7±0.5) X 


10~ 6 


E linc (keV) 


6.97 (fixed) 




ffline (keV) 


0.0 (fixed) 




Nline 


(2.4 ±0.5) x 


io- 6 


P.L. index, T 


2.36 ±0.01 




P.L. norm 


(1.51 ±0.01) 


x 10~ 3 


Emissivity index 


8 4+0-1 




Spin parameter, a 


0.998-o.oo2 




Inclination, i (deg) 


66±2 




Fe abundance (Fe/solar) 


1.0 (fixed) 




Ionization parameter, £ (erg cm s _1 ) 


1 5+°- 8 




Refiionx normalization 


(8.0 ± 1.6) x 


io- 5 


xr-dof 


1.03, 154 





NOTE. — We fit the model phabs*zedge*(power-law + zgauss + zgauss + 
kdblurcgirelionx), and fix z = 0.0257 in all relevant model components. 



Spectra and lightcurves were extracted from an 800-px 
radius circular region, with the background extracted from 
nearby source-free region of the same size to the NE of the 
source for the pn and to the SE for the MOS detectors. The 
responses were generated with rmfgen and arfgen, and the 
spectra grouped to a minimum of 25 counts per bin and a min- 
imum bin width of 1/3 of the FWHM spectral resolution at a 
given energy using the SAS 'specgroup' tool. 

2.1. Spectral analysis 

We fit the pn spectrum in XSPEC vl2, in order to deter- 
mine the energies where the power-law and reflected compo- 
nents are most dominant. Such a procedure allows us to max- 
imize the sensitivity to finding a lag between these two com- 
ponents by isolating lightcurves where the relative strength 
of each component is greatest. The model consists of an ab- 
sorbed power-law plus a bl urred reflection mo del using the 
relconv convoluti on kernel (iDauser et al.ll2010l) and refiionx 
reflection model dRoss & Fabian 20051). We use an unbroken 
power-law emissivity in relconv. In addition to this, the data 
also require an absorption edge at 0.86 keV (rest energy), con- 
sistent with an OVIII edge and two narrow emission l ines at 
approximately 6.5 and 6.97 keV (see lPorquet et alJ2007l for a 
detailed discussion of the Fe line complex in this object). We 
show the best-fitting model for t he pn spectrum in Figure Q] 
and see figure 7 of|Pprquet et al. (2007) for a ratio to a simple 
power-law. The spectrum displays a strong soft excess, which 
can be well-fit by the blurred reflection model (see Table[T]for 
spectral fitting parameters). The best-fitting spin parameter 
indicates a maximally spinning black hole. While the broad 
Fe line in this object is not especially constraining on the spin, 
it is the smooth soft excess that requires a maximally spinning 
black hole. 

Based on the shape of the spectrum, we choose the 0.3 - 
0.9 keV energy range for the soft band lightcurve, where the 
reflection component is strongest (relative to the power-law), 
and the 1.5 — 4.5 keV range for the hard band lightcurve, 
where the spectrum is dominated by the power-law compo- 
nent. Lightcurves in these energy bands were extracted using 
evselect and epiclccorr, and 20s binning. 



2.2. Frequency-dependent lags 

We first look for lags between the soft and hard bands 
as a function of Fourier-frequency in a similar manner as 
has been done for all the other sources (see references in 
secti on 1). The lightc urves exhibit significant variability 
(see iPorquet et al.l 120071 and their figure 1). The lags be- 
tween the soft and hard bands are determined from the cross- 
spectrum following the standard Fourier technique detailed in 
INowak et al.l (fT999l) . Briefly, we calculate the cross-spectrum 
for each pair of soft and hard lightcurves, averaging the cross- 
spectrum in frequency bins and also averaging the cross- 
spectra from multiple detectors. At a given frequency, /, 
the argument of the cross-spectrum gives the phase differ- 
ence between the Fourier transforms of the two lightcurves. 
This phase difference can then be converted to a time lag by 
dividing by 2nf. We therefore determine time lags in each 
frequency bin by dividing the average phase in the bin by the 
middle-frequency of the logarithmic bin. We dete r mine errors 
in the lags following equation 16 of INowak et al.l ( 1999). 

The lags between the soft and hard bands are shown as a 
function of Fourier frequency in Figure [2] The lag is defined 
such that a positive lag refers to the hard lightcurve lagging 
the soft (hence 'hard' lag), and a negative (or 'soft' lag) im- 
plies the soft lightcurve lagging the hard. We find a lag of 
-325 ±89 s in the frequency range (2-3) x 10" 4 Hz. The 
evolution of the lags with frequency shows a positive (hard) 
lag at the lowest frequencies (< 1.5 x 10 -4 Hz) transitioning 
to a negative (soft) lag at higher frequencies, (2-3) x 10~ 4 
Hz, and then becoming a zero lag at the highest frequen- 
cies (>4x 10 -4 Hz). This evolution follows the same gen- 
eral shape as has been seen in the other sources where a soft 
lag has been detected. One bin at approximately 4 x 10 -4 Hz 
shows a return to a positive lag. However, this bin is only 
2.2er from zero, and slightly larger frequency bins reduce this 
significance further. 

We model the frequency-dependent lags using two sim- 
ple transfer functions - a top-hat transfer function to fit the 
soft (negative) lags and a power-law transfer f unction to fit 
the hard (positive) lags, in a similar fashion to Zogh bi et alj 



Lags in ESO 113-G010 



3 



o 




3x10" 
Frequency (Hz) 

FIG. 2. — Frequency-dependent lags in ESO 113-G010 between the 0.3 
- 0.9 keV and 1.5 - 4.5 keV lightcurves. The lag is defined such that a 
positive lag means that the hard lightcurve leads the soft. The solid black line 
is the best-fitting transfer function, comprised of two components - a top- 
hat transfer function (blue dashed line) and power-law transfer function (red 
dashed line). 
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FIG. 3.— Coherence for ESO 113-G010 between the 0.3 - 0.9 keV and 
1,5 - 4.5 keV lightcurves. The lightcurves are coherent to approximately 
8 x 10^ Hz. 



( 120111) and Emman ouTopoulos et all ( 120111) . We find that this 
model readily reproduces the lags. Our best-fitting model 
is shown in Fig. |2] which has the top-hat transfer function 
extending from to 2390 seconds. The power-law trans- 
fer function is not well constrained as the slope and nor- 
malization are degenerate. However, any combination of 
those parameters gives the same parameters for the top-hat 
transfer function. In lag-frequency space the model for the 
hard lags is consistent with a power-law with slope of ap- 
proximately —1.7. This i s slightly steeper than the slope of 
-1 seen in other AGN (IVaughan et al.ll2003t iMcHardv etail 
l2004t lArevaloet all [2008). though a slope of -1 would im- 
ply a lag of approximately 100 s at 10~ 3 Hz, lying within the 
95.4% confidence interval derived from the data. 

We check the reliability of the lags by studying the coher- 
ence as a function of frequ ency, following the methodology of 
Vaugh an"& Nowal3 (fl997). The coherence is shown in Fig. [3] 
and remains high over the frequency range where the lags are 
detected, up to approximately 8 x 10~ 4 Hz, where the Poisson 
noise level in the hard band lightcurve is reached. The high 
coherence indicates that a high fraction of the lightcurves in 
each band is correlated and hence one lightcurve can be pre- 
dicted from the other. 




Energy (keV) 

FIG. 4. — (a) Energy-dependent lags for ESO 1 13— G010 for the frequency 
ranee (2-3) X l(r 4 Hz, calculated using the pn data only. The lags are arbi- 
trarily shifted so that the minimum lag is zero. It is the relative lag between 
each band that is important, (b) Unfolded covariance spectrum (pn data only) 
for ESO 113— G010 in the same frequency range ((2-3) X 10" 4 Hz) as the 
energy-dependent lags in the top panel. 



2.3. Energy-dependent lags and the covariance spectrum 

The energy-dependence of the lags at a given Fourier fre- 
quency can help identify the nature an d origin of the lags (e.g. 
lZoghbietalJl20Tlt iKaraet al.l 120 1 2alfbh . Here, we follow a 
similar procedure to Zoghbi et al. (201 1). Briefly, we extract 
lightcurves in narrow energy bands using the pn data only. 
Lags are then calculated between each narrow energy band 
and the reference band. For the reference band we use the 
lightcurve from the entire energy range (0.2 - 10 keV) exclud- 
ing the current narrow energy band (to avoid correlated Pois- 
son noise). This ensures a hig h signal-to-noise rat io in the ref- 
erence band. As discussed by Zog hbi et al.l d201 ll) the slightly 
different reference lightcurve used for each band leads to an 
insignificant systematic error in the lag. In Figure Ua) we 
show the lag spectrum for the frequency range (2-3) x 10 -4 
Hz where the soft lag is detected. Note that it is the relative 
lags between the bands that are important (the lags are arbi- 
trarily shifted so that the minimum lag seen is zero). As was 
shown in the frequency-dependent lags, it can be seen that the 
average lag in the soft band is lagging behind the hard band. It 
also shows that the lag is largest in the 0.4 - 0.8 keV range. If 
we modify the energy bands for the soft lightcurves slightly to 
0.4 - 0.8 keV this increases the magnitude of the lag observed 
in the (2 - 3) x 10" 4 Hz range slightly to -403 ± 83 s. 

It is interesting to note that the general shape of the lag- 
energy spectrum is quite similar to both 1H 0707-495 and 
IRAS 13224-3809 (see iKara etail 1201 2al for a comparison 
of these two objects), where the lags are seen to peak where 
the soft excess is strongest, and then drop down to a minimum 
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at 3 - 4 keV and increase again at the Fe K line. Obviously, we 
do not have the signal-to-noise in order to detect a significant 
increase in the lag above 4 keV (the lag measured in the 5 - 
10 keV range is 305 ±522 s), but, the general shape we see in 
ESO 1 13-G010 is the same. 

We also calcu l ate t he covariance spectrum (see 
IWilkinson & Uttlevl [2009b for ESO 1 13-G010 from th e 
pn data, following the methodology of lUttlev et al.l (1201 ll) . 
The covariance spectrum shows the strength of correlated 
variability between a given energy band and the reference 
band, and hence gives the spectral shape of the correlated 
variability. We calculate it in a given frequency range using 
the frequency-averaged cross spectrum in that range. The 
resultant covariance spectrum is proportional to the photon 
count rate, and so we rebin the pn response to match the 
energy binning of the lightcurves. Using the rebinned 
response, we can fit it in XSPEC. The covariance spectrum 
(unfolded using a power-law with index and normalization 
of 1) is shown in Fig. SJb). A simple power-law with index 
r = 2.1 ±0.2 fits it well, though given the signal-to-noise 
(the average fractional uncertainty is 35%), we cannot rule 
out more complex models. Note that the statistics in the 5 - 
10 keV range are too low to calculate the covariance, as is 
also demonstrated by the large uncertainty for the lag in this 
energy range. 

3. DISCUSSION 

We have found a lag of -325 ± 89 seconds in the variable 
AGN, ESO 113-G010 in the frequency range (2-3) x 10" 4 
Hz. Moreover, the evolution of the lag over the full fre- 
quency range of the lightcurve shows hard lags on the longest 
timescales changing to soft lags on intermediate timescales 
and zero lag on the shortest timescales. This evolutionary 
trend is the same as in the growing number of AGN where soft 
lags have been observed. It suggests common physical mech- 
anisms are present in all these AGN leading to hard and soft 
lags on different timescales. The soft lag can be interpreted as 
due to reverberation from the inner accretion disk, where the 
lag arises due to light-travel time between the source power- 
law emission and the accretion disk where the reflected com- 
ponent arises. Such a model predicts that the soft lag mag- 
nitude and frequency where it is observed should both scale 
with black hole mass, as the characteristic s ize-scale for the 
system is simply set by GM/c 2 . Recently, De Mar co et al.l 
(120121) observed such scaling with black hole mass using 15 
AGN where significant soft lags were detected. 

We can compare the soft lag observed in ESO 113- G010 
with what is expected from the De Mar co et al.l ([2012) black 
hole mass scaling. But first, we shoul d consider the black 
hole mass estimate for ESO 1 13-G010. iPorquet et all (120071) 
estimate the black hole mass in E SO 113-G010 using the 
mass-luminosity-timescale relation (McHar dv et al.l l2006) to 
be in the range (0.4- 1.0) x 10 7 M@, with the range due to 
the uncertainty in the bolometric luminosity. As an additional 
check on th e black hole mass e stimate, we use the optical 
spectrum of Pi etsch et al.l d 1998ft to get a separate estimate. 
Optical reverberation mapping has established a scaling rela- 
tionship b etween the broad-line region radius and A GN lu- 
minosity (iKaspi et al.l 120001 iBentz et all 120061 120091) which 
allows single-epoch black hole mass estimates when com- 
bined with the broad line width. While determined using the 
H/3 emission line, this method has been modified for vari- 
ous other emission lines too, f or inst an ce using the s t ronger 
Ha line dGreene & Hoi 120051 120071) . iPietsch et al.1 (1998) 



give the Ha line FWHM for ESO 113-G010 as 2000 km 
s _1 . From their optical spectrum of ESO 1 13-G010, we esti- 
mate the flux density at 5100A as fx = 1 x 10~ 15 erg cm" 2 s" 1 
A . Using z = 0.0257 and standard cosmology, we get XL\ 
(5100A ) = 7.5 x 10 42 erg s" 1 . We then use the iGreene 
d2005l) relationship between L(5100A) and Ln a (their equa- 
tion 1) along with the ir updated scaling with Ha FWHM from 
iGreene & Hoi ([2007) (their equation Al) to estimate a black 
hole mass of approximatel y 7 x 10 6 Mq, comp letely consis- 
tent with the estimate from Porque t et alj (120071) . 

The lag of -325 s corresponds to 9.4 R g /c for a mass of 
7 x 10 6 Mq, though dilution of the lag by any hard lag compo- 
nent, and the zero lag component that arises due to reflection 
and power-law components being presen t in both bands (see 
further discussion of dilution effects in IZoghbi et al.l 1201 It 
iKara et al.ll20l2aUbh . means that the intrinsic reverberation lag 
will be larger than this. Determining the dilution effects is 
non-trivial given that it depends on the relative strengths of the 
variable components, however, fitting of the lags with transfer 
functions, suggests that the intrinsic lag may be closer to 800s 
(see the blue line in Fig. [2]). Note, however, an added com- 
plication is that the observed lags are not adjusted for gravita- 
tional (Shapiro) time delays which would act in the opposite 
sense to dilution effects. Comparin g with the soft lag mag - 
nitude - black hole mass relation of iDe Marco et al.l (120121) . 
the lag lies significantly above the scaling relation and would 
be mor e consistent with a bla ck hole mass a factor of a few 
higher. IDe Marco et al.l ([2012) also present the scaling of the 
frequency where the soft lag is observed with black hole mass 
(the frequency decreases with increasing black hole mass). 
The frequency where the soft lag is observed (approximately 
2.5 x 10~ 4 Hz) is slightly lower than expected from the scal- 
ing, again consistent with a black hole a few times more mas- 
sive. However, it is important to note that scatter is expected 
in these scaling relations, especially given that one source, 
IRAS 13224-3809 has a soft lag that is seen to both increase 
in magnitude and decrease in frequency by a factor of approx- 
imately 3 in the high flux state compared to the low flux state 
dKara et al.ll2012al) . Thus, flux-dependence could add signif- 
icant scatter to these lag sca ling relations. Moreove r, several 
of the soft lags detected in De M arco et ail d2012l) occur at 
the lowest frequencies sampled by the lightcurve and the real 
minimum may well occur at even lower frequencies not sam- 
pled by the current data. This too can add scatter to the lag 
scaling relations. 

Another cause of scatter in the IDe Marco et al.l d2012l) re- 
lations is the black hole spin. Black hole spin changes the 
locati on of the inner disk radius dBardeen et al.lll972T: iThornd 
1974), thus objects with a maximally-spinning Kerr black 
hole will have shorter lags than objects with a non-spinning 
(Schwarzschild) black hole (assuming the rest of the geome- 
try is unchanged). Therefore, non-spinning black holes could 
be above the black hole mass-lag scaling relation, without re- 
quiring a more massive black hole. However, our spectral fit- 
ting of ESO 1 13-G010 imply an inner disk radius consistent 
with a maximally spinning black hole, suggesting that black 
hole spin is not the origin of the offset from the scaling rela- 
tions. 
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